function [CRLB,mat] = T1IRCramerRaoDT(DT, echoesmdf, noiseLvl)

% Computes the Cramer Rao Lower bound on the variance of the estimated parameters
% from measurements with echoesmdf as inversion times set.
% noiseLvl: default = 1; Specifies the noise standard deviation of the
% images (sigma in the rice distribution), the resulting CRLB depends on the noise level. 
%
% Created by Henk Smit, EMC, 01-2011 based on the work by Dirk Poot, University of Antwerp, 13-8-2007.

if nargin<3
    noiseLvl = 1;
end;

DT(2,:) = DT(2,:)./noiseLvl;

numtr = size(DT,2);
CRLB = zeros(size(DT,1),size(DT,1),numtr);
mat = zeros(3);
[A,dAdpar] = T1IRMagnitude(DT,echoesmdf);
expectIA = rice_information_ExpectVal(A);

for k=1:numtr
    FIM = (dAdpar(:,:,k)'*(dAdpar(:,:,k).*expectIA(:,k*ones(1,3)))); %HENK Fisher Information Matrix
    ws = warning('off');
    mat = inv(FIM);
    warning(ws);
    mat(:,2) = mat(:,2)*noiseLvl;
    mat(2,:) = mat(2,:)*noiseLvl;
    CRLB(:,:,k) = mat;
end;
